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Abstract. Cool giant stars arc highly dynamical objects, and complex micro-physical 
processes play an important role in their extended atmospheres and winds. The in- 
terpretation of observations, and in particular of high-resolution IR spectra, requires 
realistic self-consistent model atmospheres. Current dynamical models include rather 
detailed micro-physics, and the resulting synthetic spectra compare reasonably well 
with observations. A transition from qualitative to quantitative modelling is taking 
place at present. We give an overview of existing dynamical model atmospheres for 
AGB stars, discussing recent advances and current trends in modelling. When compar- 
ing synthetic spectra and other observable properties resulting from dynamical models 
with observations we focus on the near- and mid-IR wavelength range. 

1 Introduction 

Cool giant stars are prime targets for high-resolution infrared spectroscopy: these 
stars are luminous, emit most of their radiation at near- and mid-infrared wave- 
lengths, and they have highly dynamical atmospheres due to convection and 
pulsation, leading to strong spectral variations. The coolest and most luminous 
of these objects also have slow, massive stellar winds. They represent an impor- 
tant but not so well-understood stage in the evolution of low- and intermediate 
mass stars, and they contribute significantly to the chemical enrichment of the 
interstellar medium. These are good reasons for studying cool giants but the 
interpretation of the observed spectra requires detailed and consistent model 
atmospheres for this type of stars. 

Several of the very reasons that make pulsating asymptotic giant branch 
(AGB) stars interesting targets for IR spectroscopy, however, cause major chal- 
lenges for realistic, self-consistent modelling of their dynamical atmospheres and 
circumstellar envelopes. Many interacting processes have to be taken into ac- 
count simultaneously: propagating shock waves caused by stellar pulsation mod- 
ify the structures of the atmospheres on local and global scales, leading to strong 
deviations from hydrostatic stratification. The radiative fields are dominated by 
molecular opacities or even by dust grains forming in the cool outer layers of the 
atmospheres. Important micro-physical processes like gas-phase chemistry and 
dust formation may be severely out of equilibrium. 

Pioneering models for dynamic atmospheres of AGB stars date back more 
than two decades, but these models were severely restricted by the lack of suf- 
ficiently powerful computers and micro-physical input data. Only recently it 
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has become possible to calculate consistent dynamic models which produce rea- 
sonably realistic synthetic spectra for a wide range of stellar parameters. This 
progress is of particular interest in the context of current instrumental develop- 
ments in IR spectroscopy. 

2 Dynamic model atmospheres 

We begin this section with a definition of the term 'dynamic model atmospheres', 

in contrast to classical hydrostatic model atmospheres, stationary wind models, 
or pure circumstellar envelope models (which do not include the stellar atmo- 
sphere and wind acceleration region). The computational domain of such dy- 
namic model atmospheres is a spherical shell, extending from an inner boundary 
below the stellar photosphere to an outer boundary which can be anywhere up 
to several tens of stellar radii above the photosphere, depending on whether the 
models have a stellar wind or not. The stellar pulsation is usually simulated by 
a variable inner boundary ('piston models') and outflow of matter is permitted 
at the outer boundary in the case of a stellar wind. Within this radial range, the 
time-dependent equations of hydrodynamics arc solved, accounting for forces due 
to gravity, gas pressure and radiative acceleration of molecules and dust grains. 
Simultaneously, a more or less sophisticated treatment of radiative transfer is 
used to determine energy exchange.^ The resiilts of the numerical simulations 
are consistent spatial structures of velocities, densities, temperatures and other 
physical quantities as a function of time. Furthermore, these models may give 
mass loss rates, dust-to-gas ratios and other stellar wind properties. 

At present, four major 'strains' of dynamical models are in use. They share 
the basic ingredients given above but differ in the treatment of certain physical 
processes and the numerical methods. Listing them by their place of origin, their 
special features (and main references) can be summarised as follows: 

• The Australian models (Wood [35], Bessel et al. [3], Hofmann et al. [12]) are 
strongly linked to pulsation models, featuring the most consistent description 
of the variable inner boundary available at present (cf. Sect. 4). They have 
been used to derive synthetic monochromatic radii and molecular line profiles 
to investigate the pulsation mode of Miras (cf. Sect. 3). 

• The Iowa models (Bowcn [4]; see also Willson [30]) explore the effects of 
non-LTE cooling on mass loss, describing how thermal winds can be driven 
due to a 'calorisphere', i.e. a region heated by shock waves and kept warm by 
inefficient radiative cooling. Alternatively, in the case of rapid cooling these 
models use a parameterised dust opacity to produce outflows. 

• The Berlin models (Fleischer et al.[5], Winters et al. [31], Jeong et al. [18]) 
introduce a detailed time-dependent description of dust formation which 
results in a non-linear interaction between gas, dust, and radiation, leading 
to complex structures and variability. The models are mainly applied to 

^ A general description of the relevant physical background can be found, e.g., in 
Gustafsson & Hofner ([9]). 
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determine mass loss rates and other wind properties, especially for cool, 
dust-cnshroudcd stars. 
• The Vienna models (Hofner & Dorfi [14], Hofner et al. [16], [15]) include 
a similar detailed description of dust formation for C-rich chemistry, and 
they combine time-dependent hydrodynamics with frequency-dependent ra- 
diative transfer. The latter allows to take the effects of molecular line blan- 
keting properly into account, leading to more realistic atmospheric structures 
as demonstrated by comparison of synthetic spectra with observations (cf. 
Sect. 3). 

An overview of dynamical models, including a detailed comparison of input 
physics, is given in the recent review by Woitke ([33]). The application of wind 
models to stellar evolution is discussed by Willson ([30]). 

3 Infrared spectra: models meet observations 

When modelling the infrared spectra of cool giants, there are two basic philoso- 
phies. One approach, usually driven by new observational data, is semi-empirical 
modelling of specific observable properties of given stars. Examples are semi- 
empirical models for molecular line profiles which allow to deduce information 
about the velocity structure of atmospheres and winds (e.g. Kcady et al. [22], 
Keady & Ridgway [23]) or 'composite' models of the stellar photosphere and 
circumstellar material for the interpretation of ISO spectra, e.g. 'slab' models of 
levitated warm molecular gas (e.g. Yamamura et al. [36], Matsuura et al. [26]). 

The other method, often motivated by theoretical interests, is the construc- 
tion of more generally appliciabk^ self-consistent models, and iterative improve- 
ment of these models by comparison with observations. The models discussed in 
Sect. 2 fall into this second category: they aim at describing the time-dependent 
structure of the atmosphere of AGB stars - and in many cases the mass loss 
through a stellar wind - from first principles. In a second step, snapshots of 
these structures are used as input for detailed radiative transfer calculations, 
yielding synthetic spectra and other observable properties (colours, monochro- 
matic radii, etc.) which can be compared to observations. 

3.1 Consistency, all the way ? 

Considering the issue of consistency, it is, in principle, not sufficient if the struc- 
ture of the model is calculated in a self-consistent way. The a posteriori compu- 
tation of observable properties has to be consistent with the assumptions and 
micro-physical data used in the dynamical model. In the ideal case (which is of- 
ten realized in classical static atmospheres) the calculation of synthetic spectra 
should just be a re-computation of the radiative transfer as in the original model 
but with a higher spectral resolution. 

In practice, however, there are different degrees of inconsistency between 
physical assumptions made when calculating model structures or synthetic spec- 
tra, and it depends on the problem under investigation what is acceptable and 
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what is not. For example, when studying specific spectral lines which would not 

make a noticeable contribution to the total opacity, and therefore not influence 
the global energy and momentum balance of the model, it may be sufficient to 
introduce these lines only in the detailed radiative transfer. On the other hand, 
strong opacity sources in wavelength ranges with high flux (e.g. IR bands of 
abundant molecules like water vapour) may influence the structure of the model 
considerably, and it is essential to use the same opacity data both in the dy- 
namical calculation and in the computation of the observable properties (see 
contribution by Aringer et al., this volume). 

An even more problematic issue is the use of grey opacities in the calculation 
of the model structures, as is often done to save computing time in dynamical 
models. This may lead to unrealistic density-temperature structures, a problem 
that has sometimes been circumvented by only using the density and velocity 
from the dynamical model, and then recalculating the temperature for the com- 
putation of the spectra (e.g. Bessel et al. [3]). 

Neglecting the effects of velocity fields on spectral lines when evaluating the 
opacities during dynamical calculations may be a borderline case of acceptable 
inconsistency: in contrast to, e.g., line-driven winds of hot stars, the Doppler 
shifts are generally small in cool giants, and the molecular line blanketing prob- 
ably makes the effect even less pronounced. Other sources of uncertainty in 
this context are non-LTE effects (due to shocks), chemical non-equilibrium, or 
convective motions, which are not included in standard dynamical model atmo- 
spheres. At present, it is difficult to estimate which influence these phenomena 
have on the structures of cool atmospheres and winds. 

3.2 Crucial tests for dynamical models 

Before dynamical model atmospheres can be applied as a tool for the interpre- 
tation of observations, they have to be subjected to a number of tests which 
help to evaluate their consistency and reliability. The most simple check is, of 
course, a comparison of the hydrostatic limit case with classical model atmo- 
spheres. Testing of the full dynamical models, however, has to be done against 
other dynamical models and observations. In the latter case, it is important to 
pick a number of criteria which check different aspects of the models: 

• The global energy distribution is usually a good indicator for the overall 
structure since different parts of the spectrum are generated in different 
layers. Therefore, photometric colours and low-resolution spectra which cover 
a wide wavelength range are important for gauging the models. 

• High-resolution spectra test both the structure and dynamics of the extended 
atmosphere and wind acceleration region. They make a quantitative deter- 
mination of velocities possible and allow to distinguish between absorption 
and emission components in crowded molecular bands. 

• Time series of observations resulting from monitoring of targets give insights 
into time-dependent processes and the global dynamics which cannot be 
extracted from individual spectra. They are both the most crucial test and 
the most promising application for dynamical models. 
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Another important aspect is that these various criteria have to be fulfilled simul- 
taneously, i.e. a sequence of snapshots from one dynamical model has to fit all 
the available data for a given star, in order to claim that the model is consistent 
and realistic. 

The models currently found in the literature usually only match part of the 
criteria given above, indicating that they are probably still too specialised on 
certain aspects and do not take all necessary processes properly into considera- 
tion. In the following sections we will look at examples, focusing on the infrared 
part of the spectrum. 

3.3 Low-resolution spectra - the global picture 

The spectrometers of the Infrared Space Observatory (ISO) opened possibilities 
to obtain a simultaneous coverage of a wide wavelength range, including spectral 
regions not accessible from the ground. This provided a unique way of studying 
the global structures of atmospheres and inner wind regions of AGB stars and 
lead to a renewed interest in the mod(^lling of infrared spectra of these stars. 
After an initial phase dominated by semi-empirical modelling (e.g. Yaniamura et 
al. [36] , Matsuura et al. [26]), studies based on reasonably realistic, self-consistent 
dynamical models are beginning to emerge, both for the interpretation of ISO- 
SWS data, and ground-based observations. 

Tej et al. ([29]) discuss synthetic optical and NIR spectra resulting form 
the dynamical models of Bessell et al. ([3]) and Hofmann et al. ([12]). They 
find reasonable agreement of the overall energy distribution between 0.5 and 
2.5 /um with spectra obtained by Langon & Wood [24]. They do not, however, 
attempt to identify a single 'best fit' dynamical model for each star, since the 
sample of dynamical models is small and various spectra from different models 
may reproduce the observed spectra equally well. Furthermore, they discuss the 
problem of fitting the region around 1 fim which can possibly be attributed to 
missing or incorrect molecular opacities. 

Hron et al. ([17]) presented a first comparison of ISO-SWS spectra of R 
Scl with synthetic spectra based on the dynamical models of Hofner et al. ([16]). 
While these (grey) models were still too crude to allow for a quantitative fit, they 
could reproduce the variation of molecular features qualitatively. Since then, the 
dynamical models have been improved by including non-grey radiative transfer, 
which turned out to be crucial for obtaining more realistic structures (Hofner 
[13], Hofner et al. [15]). The observable properties resulting from this new gen- 
eration of dynamical models show a dramatic improvement compared to earlier 
grey models (see, e.g., Aringer et al. [2], Andersen et al. [1]). 

Recently, Gautschy-Loidl et al. ([8]) have performed a systematic comparison 
of spectra resulting from the non-grey models of Hofner et al. ([15]) with obser- 
vations of the C-rich AGB stars TX Psc, WZ Gas, V460 Gyg, T Lyr and S Gep 
(ISO spectra, spectra from KAO and ground-based observations, narrow-band 
colours between 0.5 and 1.1 /um). In the wavelength range between 0.5 and 5/Lim, 
they find good agreement between observations at different phases and a single 
model for each star. Figure 1 shows a comparison of observational data for TX 



6 



S. Hofner et al. 




1.0 2. 5. 

wavelength [/xm] 




wavelength [fim] 

Fig. 1. Synthetic spectra (grey) and observational data of TX Psc (black; upper panel: 
data from Joyce [20] , Lazaro et al. [25] and two KAO spectra; lower panel: ISO SWS06 
and ISO SWSOl spectra); from Gautschy-Loidl et al. ([8]). 
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Psc from various sources with synthetic spectra based on a single dynamical 
model. The 1-5 fim range (upper panel) is well reproduced but there is an inter- 
esting discrepancy between 10 and 20 fim (lower panel): the synthetic spectra 
of the pulsating but dust- and windless model show a pronounced absorption 
feature around 14 /im which is not seen in the observations. Similar problems 
have been encountered previously when fitting hydrostatic model atmospheres 
to observed spectra of C-rich AGB stars (e.g. J0rgensen et al. [19]). 

While the opposite (i.e. a missing feature in the synthetic spectrum relative to 
the observed) is easily explained with incomplete opacity data, the appearance of 
an additional absorption feature in the synthetic spectra is more astonishing, and 
several attempts have been made to explain this phenomenon (cf. J0rgensen et 
al. [19]). As discussed by Gautschy-Loidl et al. ([8]), however, dynamical models 
which show mass loss by stellar winds offer a first self-consistent explanation: 
Figure 2 compares observed spectra of the Mira variable S Cep with a dynamical 
model that has a noticeable mass loss by a dust-driven stellar wind. In this 
case, the synthetic low-resolution spectra show no absorption feature around 
14 fj,m, since the photospheric absorption feature is filled in by emission from 
circumstellar material. 

This problem of the 'missing' 14 /im feature in C-rich AGB stars illustrates 
nicely how high-resolution spectra can help to clarify a fundamental question: 
is the feature really absent (which means that there is a major problem with 
opacities and the model structures, cf. J0rgensen et al. [19]), or is it filled in by 
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emission from gas layers above the stellar photosphere, as suggested by dynam- 
ical models ? 

3.4 High-resolution spectra - atmospheric dynamics 

High-resolution spectra of AGB stars showing the variation of molecular line 
profiles with phase are an important tool for studying the dynamics of the pul- 
sating atmospheres and wind acceleration regions (cf. Lcbzcltcr, this volume, for 
an overview). In particular, vibration-rotation lines of CO in the NIR have been 
used intensively for temporal monitoring of long-period variables (e.g. Hinkle 
et al. [10]). The CO molecule has several advantages over other species: it is 
abundant, exists both in M- and C-type stars, forms at relatively high temper- 
atures deep in the atmosphere, and stays chemically inert over a wide radial 
range. Lines corresponding to different transitions arc formed in different layers, 
allowing to probe the velocity structure from the photosphere out into the wind 
region with just one molecule. 

The interpretation of the observations, however, and in particular the deriva- 
tion of absolute velocities is a non-trivial process which requires the computation 
of synthetic lines from model structures, and comparison with the observations. 
Early attempts using semi-empirical models have been mentioned above (Kcady 
et al. [22], Keady & Ridgway [23]). Here, we discuss a few recent examples based 
on self-consistent dynamical models. 

Winters et al. ([32]) studied synthetic CO fundamental and first overtone 
lines based on dust-driven wind models and compared these results to obser- 
vations of the extreme C-rich Mira IRC-l-10216. The CO fundamental lines at 
4.6 /im arc formed in the zone where the wind has more or less reached its final 
outflow velocity, showing strong P Cygni profiles and little temporal variation. 
In contrast, the first overtone lines at 2.3 /um probe the dust formation and wind 
acceleration region, and the authors interpret the observed variations in terms 
of shocks and the formation of new dust shells. They discuss the problem of 
getting a consistent fit for the strengths and shapes of the lines which indicates 
that the mass loss rate of the model is too high by about a factor of 3 while the 
dust distribution seems consistent with observations. 

Scholz & Wood ([28]) calculated synthetic CO first and second overtone lines 
and OH first overtone lines based on the dynamical model atmospheres by Besscll 
et al. ([3]) and Hofmann et al. ([12]). They derived conversion factors which relate 
velocities obtained from the Doppler profiles to the actual radial velocities di- 
rectly above and below the shock waves in the models. Applying these correction 
factors to observed velocities they conclude that Mira variables are fundamental 
mode pulsators. 

Nowotny et al. (in prep.) computed synthetic line profiles of CO fundamen- 
tal, first and second overtone lines, as well as CN 0-2 red bands (at 2/xm) using 
dynamical models of Hofrier et al. ([15]). The preliminary results show qualita- 
tive agreement with observations, e.g., reproducing the S-shaped radial velocity 
curves observed for CO second overtone lines (cf. Nowotny et al., this volume) 
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and CN, including line doubling. The amplitudes of the variations obtained with 
the models, however, are smaller than the observed variations. 

4 Current projects and trends 

In the previous seetion we have emphasised the importance of self-consistency 
when modelling the structures and spectra of cool dynamical atmospheres. We 
have mentioned a few problems of current models, such as that certain micro- 
physical processes which may be relevant for obtaining a complete and realistic 
picture are neglected for computational reasons. Apart from these 'intrinsic' 
problems, there are also a number of open issues concerning parameters and 
boundary conditions, as well as the overall geometry of the models. 

In contrast to classical stellar atmospheres which are determined by three 
stellar parameters (e.g. mass, luminosity, effective temperature) and the abun- 
dances of the chemical elements, dynamical models usually have additional pa- 
rameters describing stellar pulsation in terms of a variable inner boundary con- 
dition (see Sect. 2). There is a lack of suitable pulsation models which could 
give relations between the basic stellar parameters and quantities describing the 
pulsation, and sometimes observed period-luminosity relations are used instead. 
The ultimate goal, however, should be models that describe both the pulsation 
and the atmosphere and wind in a consistent way. A first step in this direc- 
tion has been taken by Bessell et al. ([3]) and Hofmann et al. ([12]) who used 
pulsation models to deduce boundary conditions for atmospheric models. The 
model atmospheres are a 're-computation' with higher spatial resolution and 
more complex micro-physics of the outermost layers of the pulsation models. 

A common assumption in all models discussed so far is spherical symmetry. 
This assiimption needs to be tested in the light of interferometric observations 
which seem to indicate non-spherical stellar shapes or giant spots in the photo- 
spheres of AGB stars (e.g. Karovskaet al. [21], Hofmann et al. [11]; see Scholz [27] 
for a recent review on interferometry) . Woitke et al. ([34]) investigated a possi- 
ble instability of dust formation in a non-homogeneous stellar atmosphere. This 
instability is caused by shadows which are cast by regions with a higher degree 
of condensation (and therefore higher optical depth), improving the conditions 
for dust formation in the regions which are shielded from the stellar radiation. 

Large-scale convective motions are likely to cause inhomogeneities in the 
atmospheres of cool giants. Frcytag succeeded in applying three-dimensional 
radiation- hydrodynamical 'star-in- a-box' simulations to supergiants (Freytag [6]) 
and to AGB stars, demonstrating the presence of giant convections cells (Prey- 
tag & Hofner [7]). Recently, Freytag & Hofner included a description of time- 
dependent dust formation into the CO^BOLD code to test the influence of these 
giant convection cells on grain formation. The formation and growth of dust 
grains in atmospheres of AGB stars proceeds far from equilibrium, with temper- 
ature acting as a threshold and densities of condensible material determining the 
efficiency of grain growth. Therefore, one may expect that the inhomogeneities 
in the density and temperature which are caused by the convection cells will be 
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Fig. 3. Cut through the center of a 3D radiation-hydrodynamical model ('star in a 
box') at three instants; left column: gas density; right column: dust-to-gas ratio (dark 
colours indicate low values, light colours high values; the 'stellar surface' corresponds 
approximately to the strong density gradient contour about halfway from the center 
of the box, just above the density inversion). The calculation was performed using 
171x171x171 grid points. 
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imprinted on the dust distribution around the star. Prehminary results seem to 
confirm this expectation: Figure 3 shows a cut through the center of a 3D model 
at three instants (left column: gas density; right column: dust-to-gas ratio) . The 
top panels show the model shortly after the equations describing dust formation 
have been switched on. The condensation process is fastest in the dense, cool 
gas in the lower left corner in the wake of a shock. Gradually, dust grains form 
and grow in somewhat less dense regions (middle panels) and in the wake of a 
new shock wave (upper left corner of the bottom panels). Note the rather sharp 
inner edge of the dust distribution due to the high temperatures close to the 
star. 

5 Conclusions 

Since the beginnings of time-dependent dynamical modelling more than two 
decades ago, dynamic model atmospheres for cool giant stars have overcome 
many physical and computational obstacles. Within the last few years a transi- 
tion from qualitative to consistent quantitative modelling is taking place. Cur- 
rent models include detailed micro-physics and some non-equilibrium processes 
in the computations of atmospheric structures and stellar winds. The resulting 
synthetic spectra compare reasonably well with observations, but certain con- 
sistency issues remain unsolved so far. This has to be kept in mind when using 
dynamical models for the interpretation of observations. 

In addition to the spherical models with detailed input physics, first proto- 
types of global 3D radiation-hydrodynamical models ('star-in-a-box') have been 
computed recently, investigating the dynamics of giant convection cells and their 
influence on atmospheric structure and dust formation. These 3D models are 
based on simpler micro-physics and have a lower spatial resolutions than the 
spherical models. This trade-off is necessary to keep computation times within 
acceptable limits. 

In summary, the present status of dynamic model atmospheres and the de- 
velopments that can be expected within the next few years look promising, 
regarding the interpretation of high-resolution IR spectra. 
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